Current-induced spin torques in III-V ferromagnetic semiconductors 
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We formulate a theory of current-induced spin torques in inhomogeneous III-V ferromagnetic semi- 
conductors. The carrier spin-3/2 and large spin-orbit interaction, leading to spin non-conservation, 
introduce significant conceptual differences from spin torques in ferromagnetic metals. We deter- 
mine the spin density in an electric field in the weak momentum scattering regime, demonstrating 
that the torque on the magnetization is intimately related to spin precession under the action of 
both the spin-orbit interaction and the exchange field characteristic of ferromagnetism. The spin 
polarization excited by the electric field is smaller than in ferromagnetic metals and, due to lack of 
angular momentum conservation, cannot be expressed in a simple closed vectorial form. Remark- 
ably, scalar and spin-dependent scattering do not affect the result. We use our results to estimate 
the velocity of current- driven domain walls. 
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I. INTRODUCTION 

Over a decade ago, Slonczewski and Berger— pre- 
dicted that an electrical current induces a torque on 
the magnetization of a ferromagnetic metal, and sub- 
sequent research has since identified distinct contribu- 
tions called the reactive spin transfer torque and the 
dissipative spin transfer torque, sometimes referred to 
as adiabatic and non-adiabatic, respectively* 3 . Progress 
on the understanding of this effect in metals has been 
steady ) 4 ' 5 ' 6 ! 7 ! 8 ! 9 ' 10 ' 11 ' 12 ' 13 ' 14 and the field has been stimu- 
lated by applications in spintronics and nanotechnology, 
as a way to manipulate magnetization and thus informa- 
tion. 

Current-induced spin torques arise from a small mis- 
match between the spin polarization of conduction elec- 
trons and the magnetization present throughout a ma- 
terial, and reflect the nonlocal nature of magnetization 
dynamics in inhomogeneous systems^ They are the 
converse of processes such as giant magnetoresistance 
(GMR) . The calculation of spin torques is equivalent to 
finding the conduction-electron spin density in an electric 
field, and in metals they can be easily expressed as the 
divergence of the spin current ^ Recently spin torque- 
related effects have been investigated in metal-based sys- 
tems and nanomagnets, e.g., spin-torque driven ferro- 
magnetic resonance,^ 6 spin torques in nanomagnets 17 
and in continuously variable magnetizations^ the s-d 
interaction in inhomogeneous ferromagnetsj^ 9 - current- 
induced magnetization dynamics^ and domain wall 
motion^ 

Despite these efforts, spin torques have been little stud- 
ied in materials displaying some of the richest physics and 
technological promise, namely ferromagnetic semicon- 
ductors. The most intensely studied ferromagnetic semi- 
conductors are Mn-doped III-V compounds, which fall 
into two classes: zincblende lattices, such as (Ga,Mn)As 
and wurtzite lattices, such as (Ga,Mn)N. The former 
have been much investigate d 22 ' 23 and will be the sub- 



ject of this work, while the latter are awaiting further 
experimental developments.— Recent experiments have 
succeeded in fabricating (Ga,Mn)As spin transfer torque 
devices^ while domain wall motio n 26 ' 27 and resistance^ 8 - 
have been addressed, and spin transfer physics was the 
subject of a recent review^ 9 - 

Ferromagnetic semiconductors differ profoundly from 
ferromagnetic metals. Being p-type, the carriers are 
holes described by an effective spin-3/2 and subject to 
a strong spin-orbit interaction. The hole spin is thus 
not conserved. The band Hamiltonian is the Luttinger 
Hamiltonianj^ . which, in the spherical approximation, 
is composed of a series of multipoles in spin space. 31 
The hole spin is a dipole, while the spin-orbit interac- 
tion is a quadrupole, and the precession of the spin in a 
quadrupole field is highly nontrivial. 32 Magnetism is due 
to localized Mn moments, with the exchange interaction 
mediated by itinerant holes. In the mean-field approxi- 
mation this interaction is taken into account as an effec- 
tive field (a dipole) acting on the carriers and produces a 
splitting comparable in magnitude to the Fermi energy^ 
The combined effect of this dipole exchange field and the 
quadrupole spin-orbit interaction is not a simple additive 
problem^ These materials are often in the weak momen- 
tum scattering limit (22 i.e., spT/h^> 1, where Sf is the 
Fermi energy and r is a characteristic scattering time. 
Due to the fast spin precession as a result of the spin-orbit 
interaction the relaxation time approximation is inappro- 
priate to describe spin dynamics. The carrier spin is not 
conserved which introduces arbitrariness in the defini- 
tion of the spin current. 33 Therefore, spin torques can- 
not be identified, as in metallic systems^ 4 - with the di- 
vergence of the spin current. Finally electrically-induced 
spin densities^ 4 - in bulk nonmagnetic zincblende semicon- 
ductors are forbidden by symmetry, raising the question 
of what form the spin density may take in such a system 
when the symmetry is lowered by the magnetization. 

In this paper we present a general theory of current- 
induced spin torques in zincblende III-V ferromagnetic 
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semiconductors, which, to our knowledge, is the first 
of its kind. Our study is based on the Luttinger 
Hamiltonian^S. and it assumes Eft/H 3> 1, a time- 
independent magnetization with small spatial gradients, 
zero temperature, no compensation and short-ranged im- 
purity potentials (justified by the large carrier densities 
in ferromagnetic semiconductors and the short-range na- 
ture of the exchange interaction^ 2 -). We take the spin- 
orbit interaction to be the dominant term and treat the 
exchange field in first-order perturbation theory. This is 
a good approximation for the lower end of Mn concentra- 
tions (ss 2 — 5%). We neglect also quantum interference 
effects such as weak localization. 

In recent years a number of microscopic theories 
of spin torques in ferromagnetic metals have been 
developed ! 11 ' 12 ' 13 ' 14 A common strategy is to begin with a 
uniformly magnetized state and consider small perturba- 
tions around it. The method used in this work is slightly 
different, although the physics is the same and no differ- 
ence is expected in the final results. We assume a mag- 
netization that is a function of position and has small 
gradients, and we work to first order in the gradient of 
the magnetization. We would like to note that a recent 
study has been published which considers spin torques in 
metals starting from a spin continuity equation^ 

The outline of this paper is as follows. The next sec- 
tion contains a derivation of the kinetic equation that 
will be used in the remainder of the paper. We begin 
from the quantum Liouville equation and project it in 
momentum space, then introduce so-called Wigner coor- 
dinates and derive the equation satisfied by the Wigner 
distribution. In Section III this equation is solved in the 
presence of an electric field, and in section IV the spin 
density induced by the electric field is found, which gives 
the spin torque acting on the magnetization. The form 
and implications of the results are discussed, and their 
applicability to GaMnAs is demonstrated. Finally, the 
domain wall velocity as a result of the spin torque is es- 
timated in Sec. |V] 



II. KINETIC EQUATION 

The typical setup for a spin torque experiment consists 
of two slabs of ferromagnetic material with non-collinear 
magnetizations, separated by a tunnel barrier. Since the 
magnetizations of the two slabs are non-collinear there is 
a region near the interface over which the magnetization 
changes. To determine the continuum limit of this setup 
one can begin by visualizing a large number of slabs put 
together, with slight variations in the direction of the 
magnetization of each slab. Then one can imagine the 
interfaces between the slabs disappearing, leaving one 
large sample with an inhomogeneous magnetization, in 
such a way that the gradient of the magnetization varies 
little over distances comparable to the lattice spacing. 
The gradient expansion in the magnetization that fol- 
lows from this procedure is valid as long as the length 



scale on which the magnetization varies is much longer 
than the relevant length scales of the carriers, i.e., the 
Fermi wavelength and mean free path. 

Spin torques appear when an electrical current flows 
through such a material. In a ferromagnetic semiconduc- 
tor the magnetization is a result of the Mn ions which 
interact by means of the exchange coupling mediated by 
itinerant holes. The holes themselves have a spin po- 
larization, and in equilibrium the hole spin polarization 
follows the magnetization. When an electrical current 
flows through the sample one can think for example of 
a hole which is taken from position r, where its spin is 
parallel to the local magnetization at r, and transporting 
it to r + Sr. The magnetization at r + Sr is slightly dif- 
ferent from the magnetization at r, so a torque is exerted 
by the itinerant hole on the magnetization. What makes 
the situation in ferromagnetic semiconductors more diffi- 
cult and more interesting is that while the hole is moving 
from r to r + Sr it is subject to the strong spin-orbit in- 
teraction, which acts to randomize its spin. In order to 
determine the effect of spin-orbit interactions, which are 
wave vector-dependent, on the itinerant holes and ulti- 
mately on the inhomogeneous magnetization we need to 
study the kinetic equation, which takes into account both 
the momentum dependence and the position dependence. 

In this section we will derive a kinetic equation suit- 
able for describing inhomogeneous ferromagnetic semi- 
conductors in an electric field. We consider the system 
to be described by a density operator p, which obeys the 
quantum Liouville equation 

S + ilM-a (i) 

The total Hamiltonian H contains contributions due to 
the band Hamiltonian H Vl the scalar impurity potential, 
the exchange interaction between delocalized holes and 
localized Mn moments, and the electric field. These will 
be given below. The Liouville equation is projected onto 
a set of states \v,k s ) of definite wave vector k and spin s, 
which are assumed to be Bloch functions and eigenstates 
of the Luttinger Hamiltonian H v . The matrix elements of 
p in this basis are p ss '(k,k') = p(k,k') and are treated 
as matrices in spin space. H v is diagonal in fc and its 
matrix elements in this basis are H v = iJ„(fc) 30 i 31 

^ = |l[7ifc 2 +7(!fc 2 -2(fc-S) 2 )], (2) 
2m 

where S is a vector of spin-3/2 matrices. The term pro- 
portional to the Luttinger parameter 71 gives the hole 
kinetic energy. For k 7^ 0, the term proportional to 7 
separates the heavy hole (HH) and light hole (LH) states, 
i.e., it is the spin-orbit coupling that plays a central role 
in the present analysis. The Mn 2+ ions give rise both to 
a net magnetic moment, through the hole-mediated ex- 
change interaction, and to scattering, which has a scalar 
part and a spin-dependent part. These are contained in 
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the Hamiltonian 

H M n(r) = H r + V(r - Ri) a t ■ S] , (3) 

where the sum runs over the positions Ri of the Mn 2+ 
ions, with s/ the Mn spin. We approximate the interac- 
tions represented by Hmu as short-ranged, so that U(r — 
R I ) = U6(r-Ri) and V(r-Jlj) = ( J pd /V) S(r - R T ), 
with J p d the exchange constant between the localized Mn 
moments and the itinerant holes, and V the sample vol- 
ume. The matrix elements of i?Mn in the basis { |iifcs) } 
are decomposed into a part H v d diagonal in k, which 
gives the net magnetization M, 

Hft/ = H pd = s-S^M-S, (4) 

where iVMn is the number of Mn 2+ ions and we assume 
all Mn spins polarized setting sj = s, and a part off- 
diagonal in fc, which causes spin-dependent scattering 
and will be given below. We will concentrate in this work 
on the case in which the exchange splitting is smaller 



than the spin-orbit coupling at the Fermi energy, i.e., 
\M\ <C 2jh 2 kp/m, where kp is the Fermi wave vector, 
and we work to first order in m\M\/ (2jh 2 k F ). 

In studying inhomogeneous magnetization one must 
account for real-space as well as wave- vector dependence, 
which is accomplished by defining a Wigner distribution. 
The Wigner function corresponding to the one-particle 
density matrix p(k, k') is 

A(r) = |^ e iQ >(g + Q/2,<7-Q/2). (5) 

where q = (k + k')/2 and Q = k — k'. The next step is to 
derive an equation describing the time evolution of the 
Wigner distribution f q (r). The kinetic equation for the 
Wigner function / = f q {r) is obtained by projecting the 
quantum Liouville equation onto the basis states \uk s ), 
then making the transformation (0. The Hamiltonian 
H is diagonal in wave vector. The first step, in which no 
approximations have been made, gives us the Liouville 
equation in terms of the so-called Wigner coordinates q 
and Q 



^4 I wf el *' r{Hq+Pq+q - =4? / ^ e<Qr (^-- -/w-.)- (6) 

To obtain a transparent kinetic equation it is necessary to expand H q± around the wave vector q, which requires 
some care. In this article we are working in a basis in which the functions depend on the wave vector q, and we 
require a formulation of the kinetic equation that is manifestly gauge covariant. This means that the kinetic equation 
should not acquire additional terms if the basis functions are subjected to a q-dependent rotation. The end result 
of this requirement is that ordinary q-derivatives are replaced by covariant derivatives (r-derivatives, denoted by 
V, remain unchanged since the basis does not have position dependence.) The covariant q-derivative is defined by 

= jf- — i [72., /], where the gauge connection matrix 72. SS ' = (u qs \i g^ 3 ' )- All this implies that in our derivation 

we must replace the ordinary derivatives by covariant derivatives. Expanding H q± w H q ± ^ • -j^ equation ([6]) can 
be written as 

f + It + ■ - 4 £/ 0«^(CWW- - (7) 



where {} denotes an anti-symmetrized dot product, that 
is the symmetrized scalar product between vectors a and 
b is given by {a ■ b} = a ■ b + b ■ a. The scattering term 
represented by the RHS is dealt with in the appendix. 

In a constant and uniform electric field the total Hamil- 
tonian contains an additional term containing the elec- 
tromagnetic potential V — eE f, where r is the position 
operator. This term is diagonal in real space. Following 
the spirit of the derivation presented above, this term 
appears on the right side of the kinetic equation in the 
same way as U, and is expanded as follows 



The spatial gradient of the external electrical potential is 
equal to the electric field — W = E. In this work we will 
be studying the response of the system to linear order in 
the electric field. 

When formulating a kinetic equation, which takes into 
account the variation of the Wigner function in real space 
as well as in momentum space, it is necessary to sin- 
gle out the length and wave-vector scales relevant to the 
problem under study. In the work at hand we consider 
carriers which are delocalized in real space and are de- 
scribed by Bloch states, for which the wave vector is a 
good quantum number. Nevertheless it must be borne in 
mind that the carrier occupies a finite range of real and 
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momentum space, denoted by Ar and Aq respectively, 
which are determined in such a way as to be consistent 
with the Heisenberg uncertainty principle. In the course 
of a scattering event in which a carrier with wave vec- 
tor q interacts with the potential of an impurity and its 
wave vector changes from q to k, it is necessary as well 
as physical to assume that the wave- vector spread Aq as- 



where S_e = —eE- (Df /Dq) is the covariant form of the 

usual source term due to E. The term J(f) represents 
the scattering term, which is discussed in detail in the 
appendix. The scattering term takes into account the ef- 
fect of the potential U, which represents the part of the 
Hamiltonian Hmu which is off-diagonal in wave vector. 
An explicit form for the scattering term will be given be- 
low when we discuss the solution of the kinetic equation 
in an electric field. 



III. SOLUTION OF THE KINETIC EQUATION 

The equilibrium distribution / eq is the solution to Eq. 
© in the absence of external fields, Yje = 0. To lead- 
ing order in \U\ 2 this solution is / cq = fo(H v + H p d), 
with /o the Fermi-Dirac function. It is straightforward 
to check that this form of the Wigner function satisfies 
the kinetic equation ^ when the RHS is equal to zero. 
The form of his solution shows that in equilibrium the 
spin polarization of the holes follows the magnetization 
of the Mn, which is contained in the exchange part of the 
Hamiltonian H pe i- 

Next, in the linear response regime, we search for a 
solution of the kinetic equation for nonzero Eg, which 
will yield the spin density induced by E. Since the spin 
density induced by the electric field will be a function of 
position and will in general not be parallel to the local 
magnetization, this will immediately give the spin torque 
exerted by the conduction holes on the magnetization. 
The method we use to solve the kinetic equation is as 
follows. First, we divide every matrix M. in the problem 
into M m + M out , where M m has elements only within 
the HH and LH subspaces, while M° ut has matrix ele- 
ments only between these subspaces. Schematically this 
can be summarized by 

M= (inout\ 

\out in I - ' 

One compelling advantage of this decomposition is that 
commutators and anticommutators of matrices belonging 
to either the in or out sectors do not mix these sectors. 



sociated with the carrier size is m! uch smaller than the 
typical momentum transfer in scattering processes k — q. 
Furthermore, it is assumed that the magnetization M 
varies over length scales much larger than interatomic 
separations. With these assumptions, the kinetic equa- 
tion in an electric field E takes the form (in agreement 
with the form found by Carruthers and Zachariasen^i) 



(9) 

I 

The following list covers all the possible combinations of 
commutators and anticommutators of matrices belonging 
to either the in or the out sectors 



[in, out] 


= out 


[in, in] 


— in 


[out, out] 


= out 


{in, out} 


= out 


{in, in} 


— in 


{out, out} 


= in. 



Another advantage of this decomposition is that it aids 
us in constructing a physical picture of spin torques and 
their relation to spin precession. The decomposition into 
an in and an out sector in effect singles out spin pre- 
cession as a result of the spin-orbit interaction. The in 
sector of the density matrix represents spins that are sta- 
tionary under the action of the spin-orbit interaction, or 
alternatively the fraction of the spins that are in eigen- 
states of H v . The out sector on the other hand rep- 
resents spins that precess under the action of the spin- 
orbit interaction. This decomposition determines which 
spin torques are due to the hole spin precession, which, 
unlike the precession of spin- 1/2 electrons, cannot be at- 
tributed to an effective magnetic fields Being in the 
weak momentum scattering regime EftIYl 3> 1, we do 
not consider scattering in the out sector or between the 
in and out sectors (it can be shown that both of these 
terms yield corrections linear in \U\ 2 ). The Wigner func- 
tion / has two parts, f m in the in sector and f out in the 
out sector, and the kinetic equation is broken down into 
two coupled equations for f m and f out 

Bf in i » ■ 

-^ r + j:[H m ,.r] + J(f m ) = + (12a) 

f) fOUt A 

J i r zjtn rout-\ -s^out , \^out /i ou\ 

~~dt~ h'- * = E gr ' ^ ' 

There are two source terms in each equation, namely E^ 1 
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and S*™ in the in sector, and E^ M * and E°"* in the out 
sector. To obtain these source terms one needs to ex- 
pand all quantities in the gradient of the magnetization 
and keep terms to zeroth and first order in this gradient. 
To zeroth order in the gradient of the magnetization the 
source terms are j]™/ otlt = (eE/K) ■ (-D/ cq /Z?q) m ^° u *, 



which are found simply by taking Y,e defined above and 
substituting f eq for the Wigner function. When the ex- 
pansion is continued to the next order, the source terms 
linear in the gradient (gr) of the magnetization are 



1 

2h 
1 

2fi 



VH° ut 
VH in ■ 



Df c 



Dq 

jout 

Dq 



i r DH out 

2h \ Dq 

i r DH in 



2h { Dq 



^ jout 

^ £OUt 



1 

2h 



DH l 



Dq 



(13a) 
(13b) 



To obtain Eq. ([13]) we have assumed a small spatial gradient VH = VH pd implying a small variation SM -C |(Af)|, 
and we worked, as stated, to first order in m\M\/(2 ; yh 2 kp). After some simplification we obtain for the scattering 
term acting on f m 



fin _ fin p2 1 m * 

J{D = — + — f n - - r s f™ r s + 

T T t Tq z h z 



47T 



' , 9f in 

(f m - f ' m ) (H p d - H' pd ) - q —q^- (H P d - H p 



l pd) 



, (14) 



where the bar is an average over directions in momentum 
space, t _1 = N M n\U\ 2 m*q/(Virh :i ), f = f(q,8,<p) and 
/' = f{q,0', (/>'), and 4> are the polar and azimuthal 
angles of q (analogously for q'), and m* is the carrier 
effective mass, which is to/ (71 — 27) in the HH subspace, 
and m/(7i + 27) in the LH subspace. 

For simplicity and without loss of generality we choose 
E || y and (M) || z so that M x . y {r) < M z (r). We 
solve Eqs. (TT^j) as follows: the equation for f out is first 
solved with S™* as the initial source, and the solution 
/° ut thus obtained is substituted into T,°f and E™. The 
equation for f m is solved in an analogous fashion. The 
solutions to the equations for f m and f out involve expres- 
sions of the form e i- HZ,lt / h j\/[ e ~ lHmt / h ^ an d f out is easily 
found. This is because in the out sector the product 
e iH" l t/hj^out e -iH zn t/h con t a i ns on ly functions of time of 
the form sinajt and coswi, with lu = 2ftrjq 2 /m the en- 
ergy difference between the HH and LH bands when the 
magnetization is zero. The steady-state solution for f out 
therefore involves only a straightforward time integral of 
the kind customarily encountered in linear-response the- 
ories. The equation for /*" takes more effort due to the 
presence of the scattering term and we only summarize 
the method here (it is described in detail in a recent pub- 
lication by two of us2£) . The in sector represents the part 
of the Wigner function that is stationary under the action 
of H v . Nevertheless, the full Hamiltonian is H v + H p d, 
and the commutator [H p ^,f m ] is not zero. In a manner 
similar to the decomposition of / into f m and f out , f m 
itself is split into a part that commutes with H™ d -, and 
a part that does not. It can be shown^ that the com- 
muting part yields a correction to the Wigner function 
that is linear in r while the non-commuting part gives a 



correction that does not depend on r. However, we find 
that all contributions to / average to zero over directions 
in momentum space except /g r u *. This implies that all 
contributions from f m average to zero over directions in 
momentum space. /°"* gives rise to a spin density <S 
that is independent of scattering. It is discussed in detail 
below. 



IV. SPIN TORQUES 

The only contribution to the spin density in an electric 
field comes from The three components of the spin 

density S that this correction to the Wigner function 
yields are 




These equations are the central result of our work. The 
dimensionless quantities r\i and C«j with i = x,y,z are 
functions of the Luttinger parameters 71 and 7. For 
GaMnAs we find (all xl(T 4 ) J] x = n z = 3.66, r\ y = 5.52, 
Cx = 11-56 and £ y — 6.16. The steady-state spin den- 
sity is not collinear with the magnetization, so there will 
be a torque on the magnetization giving a precession fre- 
quency of magnitude J p d\S\. Taking p = 1.2 xlO 20 cm~ 3 , 
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E y = 100 kV/m, and estimating the change in the mag- 
netization as 20% over one lattice spacing, the time scale 
of this precession is 200 ns - less than in metals, but M 
itself is also typically one order of magnitude smaller. 



A. Discussion 

The fact that the spin torque comes only from /°™* 
implies that the steady-state spin density is due to pre- 
cession under the action of both the spin-orbit interac- 
tion and the exchange field. The fraction of the spins 
that is conserved, which would yield a term oc r, gives a 
contribution that averages to zero in momentum space. 
The quantities r\% and Q decrease with increasing spin- 
orbit interaction (given by 7), suggesting the spin-orbit 
interaction reduces the spin torque. This agrees with the 
finding that there is no electrically-induced spin density 
in the corresponding nonmagnetic systems^ i.e. in the 
limit of large spin-orbit interaction 7. This limit is equiv- 
alent to restoring the spherical symmetry of the Luttinger 
Hamiltonian of Eq. (J2J) , which in ferromagnetic semicon- 
ductors is broken by the magnetization. 

An important difference from ferromagnetic metals is 
that, in Eq. (JTSJ) , there is no contribution from scatter- 
ing, either scalar or spin-dependent. This fact indicates 
that the dominant spin torque in ferromagnetic semicon- 
ductors in the weak momentum scattering limit is intrin- 
sic. This observation agrees with the results of Jungwirth 
et al.r^- who studied the anomalous Hall effect in ferro- 
magnetic semiconductors in the regime Eft/Ii 3> 1 and 
found similarly that the role of scattering is secondary. It 
is also related to the absence of electrically-induced spin 
polarization in bulk nonmagnetic zincblende materials. 
Generally, such a spin polarization is due to the fraction 
of spins that is conserved 3 -^ and is linear in r, but this 
spin polarization is forbidden by symmetry in zincblende 
lattices.^ The magnetization breaks the cubic symme- 
try of the lattice and gives a steady-state spin density, 
but the term linear in r still averages to zero. We come 
back to the comparison of our result with result found 
for ferromagnetic metals in the next section. 

We find that an electric field E \\ x corresponds to 
the permutation x <-> y in Eq. (|15p . Yet for a given ori- 
entation of E, unlike in ferromagnetic metals, in ferro- 
magnetic semiconductors there is no symmetry between 
the different components of the spin density for the fol- 
lowing reason. In metals spin is conserved and spin 
torques can be derived phenomenologically directly from 
the Landau-Lifshitz-Gilbert equation (the so-called book- 
keeping argument&ii) . One assumes an itinerant spin 
passes a localized moment at r, lines up with it, then 
moves on to another moment at r+dr and exerts a torque 
on this moment. This relates M(r + Sr) to M(r) and 
gives a simple vector-product form for ,S l (r) i 11 i 12 ' 13 i 14 In 
ferromagnetic semiconductors the spin-orbit interaction 
acts to randomize the itinerant spin moving between r 
and r + 5r, and there is no simple relationship between 



M (r + 5r) and M(r). Such a book-keeping argument 
is thus not valid and there is no symmetry in the final 
expression for the spin density. 

We would like to comment on one last aspect of the 
relationship between the hole spin polarization and the 
magnetization in ferromagnetic semiconductors. The cal- 
culation presented in this work relies on a mean-field de- 
scription of the magnetization and hole spin polariza- 
tion. In this picture the itinerant holes are subject to 
an average magnetic field due to the Mn 2+ ions, and the 
Mn 2+ ions in turn are subject to the itinerant hole spin 
polarization, which can also be regarded as an average 
magnetic fields Since it is assumed that the spin-orbit 
interaction has spherical symmetry, there is no easy axis 
for the magnetization in the absence of an electric field. 
However, once the electric field is applied it is natural to 
ask whether the direction of the electric field provides an 
easy axis for the magnetization, in other words whether 
the magnetization in the direction of the electric field in- 
creases. We find that this indeed is true, but the increase 
in the magnetization is second order in the ratio H p <i/£f 
and is not significant. 

In ferromagnetic metals, in which spin-orbit coupling is 
negligible, angular momentum is conserved. As a result 
spin torques in these materials can be encapsulated into 
a set of simple, compact, rotationally-invariant vectorial 
expressions. In ferromagnetic semiconductors, in which 
spin-orbit interactions are usually the dominant energy 
scale, angular momentum is not conserved and the final 
expressions for the spin torques cannot be expected to 
have rotational invariance. In principle spin-orbit inter- 
actions, which couple the spin and the lattice, should give 
magnetic anisotropy and anistropic spin torques as well. 
The anisotropy in our result for the spin density is thus 
a direct result of the intrinsic spin-orbit interactions. 

B. Parameters and applicability for GaMnAs 

We shall assume a doping density riMn = J) = 1.2 x 
10 20 cm -3 , corresponding to x = 2.2%, J p d — 54 meV 
nm 3 as discussed in Ref. [22j and the lattice constant 
a = 5.6533A The Fermi energy is found as 

,2m h e F , 3/2 / 2771/6^,3/2 2 



e F = 1.633- — (3tt 2 t7) 2/3 = 2.1 x 10~ 20 J 

2mo 

and the heavy and light hole Fermi wave vectors are kh = 
1.43 x 10 9 m _1 and fc ; = 0.55 x lOW 1 . The heavy hole 
and light hole masses are = 0.538 x 10 _30 kg and mi = 
0.076 x 10 _30 kg. These numbers also give the magnitude 
of the effective field \H p d\ — nun Jpd < Smu >= 2.52 x 
10~ 21 J, meaning that the ratio \H p d\/£F = 0.12, so it is 
safe to do perturbation theory. 

We also want to work out £FT p /h. The Fermi energy 
is 2.1 x 10~ 20 J, which means £ft p /H > 1 for any mo- 
mentum scattering time t p > 5 x 10~ 15 s. For example 
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for £ pTp I h w 10 we require t p = 5x 10~ 14 s, correspond- 
ing to a light-hole mobility of approximately 1000cm 2 /Vs 
and a heavy-hole mobility of approximately 200cm 2 /Vs. 
Thus the theory is on very firm ground even for extremely 
low mobilities. 



V. DOMAIN- WALL MOTION 

As an application of the central result in Eq. (| 1 5[) we 
calculate the spin torque on a domain wall, and the re- 
sulting domain-wall velocity. We choose the current and 
variation of magnetization in the y-direction. Further- 
more, we use rj x = r\ z such that Eq. (|15p reduces to 



o — 



-3/2 



dy 



M 



Tlx 



1 Myf 



(17) 



with y a unit vector in the y-direction. The spin-transfer 
torque that acts on the magnetization is given by 



dM 
~dt 



J-pd 

h 2 



MxS. 



(18) 



Using Eq. (|15p we rewrite this as an equation for a unit 
vector n in the direction of magnetization, i.e., M = 
nMnJpdSMnQ, with Smu = 5/2 the spin of one Mn atom. 
We find that 



~dt 



-vtt x 



d_ 

dy 



ft 



current 

with the velocity v given by 



1 Q y y 



(19) 



v 



nMneEm 1 / 2 r) x J 2 d SMn 

fc2 3 / 2 



The result for the current-induced torques in Eq. (JT9J) 
has the form of an anisotropic dissipative spin transfer 
torque;— The reactive spin transfer torque contribution 
is equal to zero. These results are understood by noting 
that we have considered strong spin-orbit interactions, 
and that have done perturbation theory in the magneti- 
zation. 

It is common to define a dissipative coefficien t 11 ^ 3 (3 
such that v ~ with j the current density. Because 
our result for the spin transfer torque is independent of 
r, and because j ~ r, we would find that (3 ~ 1/r, i.e., 
resistivity-like. This is somewhat surprising as recent 
studies 41 indicate that the Gilbert damping constant etc, 
which is believed to be similar though not exactly equal 
to /3, predominantly has intra-band contributions that 
are conductivity-like. However, a direct comparison is 
not possible because in the present paper we perform an 
expansion in the magnitude of the magnetization whereas 
Ref. [4l| calculates a G by determining the transverse re- 
sponse function. 



The velocity v divided by the Gilbert damping con- 
stant provides an estimate for the domain wall velocity 
so that 



X . 



(21) 



Although Sinova et al. do not explicitly consider the 
regime of parameters quoted in the previous section, their 
calculations 4 ^ (see also Ref. [4l|) suggest that the Gilbert 
damping is very small a G ~ 0.001 in this regime. Using 
this result we find that X ~ 1 m/s, in agreement with 
experimental results for the domain-wall velocity^ 

To investigate more quantitatively the effect of the 
anisotropy in the spin transfer torque, determined by the 
ratio rjy/r] x , we consider specific model for a magnetic 
domain wall. We consider a thin film, in which there is 
a constant hard-axis anisotropy K± perpendicular to the 
film and an easy-axis anisotropy Within the model 
for a rigid domain wall proposed by Tatara and Kohnc 42 
(see also Ref. [HI), the domain wall is described by two 
collective coordinates: the position X(t) and the chirality 
4>o(t). The chirality is the angle with which the magnetic 
moment in the center of the domain wall tilts out of the 
easy plane. Using the results from Ref. [42j and Ref. [43[ 
we find the equations of motion for the domain- wall col- 
lective coordinates. They are given by 



X 



K 



a G (p G 



■ sin 20o 



4>o + Q'G 



X 



A 



(l + (5cos 2 (j) ) 



(22a) 
(22b) 



where A = y J/K-g is the width of the domain wall. 
Note that S goes to zero for r\ x — > r) y . Note that, in 
addition usual dissipative spin transfer torque contribu- 
(20) tion to these equations discussed in earlier work , 42 i 43 we 
find a chirality-dependent anisotropic contribution pro- 
portional to (5. 

The above equation can be solved analytically. From 
this we obtain the average drift velocity as a function of 
the applied electric field, as shown in Fig. [1] From this 
figure, we observe a Walker-breakdown-like behavior, 44 
i.e., the domain wall velocity reaches a maximum and 
then becomes smaller. Physically, the breakdown is due 
to the transition of rigid motion to precessional motion 
of the domain wall and is well-known from field-driven 
domain-wall motion. Our results are understood from 
the fact that dissipative spin-transfer torque enters the 
equations of motion for the domain wall in the same way 
as an external magnetic field. Note that the anisotropy 
6 alters the result for the domain- wall velocity somewhat 
with respect to the isotropic (<5 = 0) situation, but plays 
no qualitatively important role. 

As a final remark, we note that in the calculations 
presented here we have neglected the effects of finite 
temperature^ 3 - and pinning of the domain wall. This, in 
addition to the fact that the experiments of Yamanouchi 
et alw^ are in a different regime of doping than consid- 
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(1 +a 2 )tiX 




the case of strong spin-orbit coupling. We have applied 
our results to the case of current-driven domain wall mo- 
tion and have estimated the resulting domain-wall ve- 
locities. We find domain-wall velocities that are of the 
same order of magnitude as experiments, although the 
available experimental results^ are in a different regime 
of parameters than considered in this paper. Therefore, 
a more quantitative comparison between theory and ex- 
periment is at present not feasible. 



FIG. 1: Average drift velocity as a function of the applied 
electric field and for 5 = 1/3 (solid line) and S — (dashed 
line) . 

ered here, makes a direct quantitative comparison not 
possible. 

VI. CONCLUSIONS AND 
ACKNOWLEDGEMENTS 

In conclusion, we have established a microscopic theory 
of spin transfer in III-V ferromagnetic semiconductors for 



We wish to acknowledge enlightening discussions 
with A. H. MacDonald, M. Tsoi, G. E. W. Bauer, 
L. W. Molenkamp, G. Schmidt, B. J. van Wees, M. van 
Veenendaal, D. J. Keavney and Zhuge Liang. The re- 
search at Argonne National Laboratory was supported 
by the US Department of Energy, Office of Science, Of- 
fice of Basic Energy Sciences, Contract No. DE-AC02- 
06CH11357. The research at Utrecht University was sup- 
ported by the Netherlands Organization for Scientific Re- 
search (NWO) and by the European Research Council 
(ERC) under the Seventh Framework Program (FP7). 



1 J. Slonczewski, J. Mag. Mag. Mat. 159, LI (1996). 

2 L. Berger, Phys. Rev. B 54, 9353 (1996). 

3 D. C. Ralph and M. D. Stile s, larXiv:071 1.46081 
(2007); P. M. Hane y et al, larXiv:0709.3862l (2007); 
Y. Tserkovnyak et al. larXiv:0709.2937l (2007). 

4 Ya. B. Bazaliy et al, Phys. Rev. B 57, 3213 (1998). 

5 M. Tsoi et al, Phys. Rev. Lett. 80, 4281 (1998). 

6 M. D. Stiles and A. Zangwill, Phys. Rev. B 66, 014407 
(2002). 

7 J.-Ph. Ansermet et al, IEEE Trans. Mag. 40, 358 (2004). 

8 Z. Li and S. Zhang, Phys. Rev. Lett. 92, 207203 (2004). 

9 A. Thiaville et al, Europhys. Lett. 69, 990 (2005). 

10 S. E. Barnes and S. Maekawa, Phys. Rev. Lett. 95, 107204 
(2005). 

11 R. A. Duine et al, Phys. Rev. B 75, 214420 (2007). 

12 F. Piechon and A. Thiaville, Phys. Rev. B 75, 174414 
(2007). 

13 H. Kohno et al, J. Phys. Soc. Jpn 75, 113706 (2006). 

14 Y. Tserkovnyak et al., Phys. Rev. B 74, 144405 (2006). 

15 Y. Tserkovnyak et al, Rev. Mod. Phys. 77, 1375 (2005). 

16 J. N. Kupferschmidt et al., Phys. Rev. B 74, 134416 (2006). 

17 X. Waintal and O. Parcollet, Phys. Rev. Lett. 94, 247206 
(2005). 

18 J. Xiao et al, Phys. Rev. B 73, 054428 (2006). 

19 A. Rebei et al, Phys. Rev. B 72, 064408 (2005). 

20 J. Fernandez-Rossier et al, Phys. Rev. B 69, 174412 
(2004). 

21 G. Tatara et al, J. Phys. Soc. Jpn. 75, 064708 (2006). 

22 T. Jungwirth et al, Rev. Mod. Phys. 78, 809 (2006). 

23 J. Konig et al. in Electronic Structure and Magnetism of 



Complex Materials, Springer Series in Material Sciences 
54, (Springer, 2003); J. Konig et al. in Interacting Elec- 
trons in Nanostructures, Lecture Notes in Physics 579 
(Springer, 2001). 

24 D. J. Keavney et al, Phys. Rev. Lett. 95, 257201 (2005). 

25 J. Wunderlich et al. , larXiv:0707.3329 1 (2007) . 

26 A. K. Nguyen et al, Phys. Rev. Lett. 98, 146602 (2007). 

27 M. Yamanouchi, D. Chiba, F. Matsukura, T. Dietl, and H. 
Ohno, Phys. Rev. Lett. 96, 096601 (2006). 

28 A. K. Nguyen et al, P hys. Rev. Lett. 9 7, 136603 (2006). 

29 H. Ohno and T. Dietl. larXiv:0712.3247l (2007). 

30 J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 

31 N. O. Lipari and A. Baldereschi, Phys. Rev. Lett. 25, 1660 
(1970). 

32 D. Culcer et al, Phys. Rev. Lett. 97, 106601 (2006). 

33 J. Shi et al, Phys. Rev. Lett. 96, 076604 (2006). 

34 E. L. Ivchenko and G. E. Pikus, JETP Lett. 27, 604 (1978). 

35 G. Tatara and P. Entel, Phys. Rev. B 78, 064429 (2008). 

36 H. Ohno and T. Dietl, J. Magn. Magn. Mat. 320 1293 
(2008). 

37 P. Carruthers and F. Zachariasen, Rev. Mod. Phys. 55, 
245 (1983); V. P. Sillin, Sov. Phys. JETP 6, 945 (1957). 

38 D. Culcer and R. Winkler, Phys. Rev. B 76, 245322 (2007). 

39 T. Jungwirth et al, Phys. Rev. Lett. 88, 207208 (2002). 

40 Jairo Sinove, et al, Phys. Rev. B 69, 085209 (2004). 

41 Ion Garate and Allan MacDonald, larXiv:0808. 3923vl 
[cond-mat .mtrl-sci] . 

42 G. Tatara and H. Kohno, Phys. Rev. Lett. 92, 086601 
(2004); 96, 189702 (2006). 

43 R.A. Duine, A.S. Nunez, and A.H. MacDonald, Phys. Rev. 



9 



Lett. 98, 056605 (2007). 
44 N.L. Schryer and L.R. Walker, J. Appl. Phys. 45, 5406 
(1974). 



APPENDIX A: SCATTERING TERM 



The scattering term J(f) is 



U,e->"< [[/,/] 



(Al) 



with rj a regularization factor and the impurity average 

( H ul H Mn)^<i' = N Mn \U\ 2 (l + T s ), where AT M „ is the 
number of Mn impurities and T s = J^ d (s • S) 2 / (\U\ 2 V 2 ). 
The derivation of this general form of the scattering term 
is discussed in a recent paper— and the notation will 
be explained in detail below. In terms of the Wigner 
distribution the scattering term can be expressed as 



TlMn f d 3 K 

n 2 J (2tt)3 j t , 

tlMn f d 3 K 



dt' eT rit ' u, 



-iH K t' , 



Unqfq f K.U nq) 



iH a t' 



(A2) 



r 



We must note that, in the approximation we are using, at the explicit form of this term and note that, because 
the scattering term acts only on f m , which brings about it involves only f m , this term commutes with the time 
some simplifications. These become apparent if we look evolution operators. 



J 



J(/ m ) 



TlMn 

TlMr, 
If 



d 3 K 

(2^)3 

d 3 n 
(X) 3 " 



dt'e^ 1 ' U(e- lH ^'Ue lH " t ' f q n - /**e -ar - t 'l7e <fl « t ') 



U. 



(A3) 



In the approximation used in this paper, the Hamiltonian entering the scattering term is the projected 2x2 Hamilto- 
nian for each subspace. The scattering potential has two parts, one a scalar and one which is spin-dependent. Taking 
into account also the exchange splitting of the bands, there are three contributions to the scattering term: scalar 
potential + kinetic energy (= Jo), spin-dependent potential + kinetic energy (= J m ), scalar potential + exchange 
energy (= J s ). The former two sum up to 



J (/ m ) + Jm(D 



1 

r 

fin 
Jq 



2-KUMr, 



n Mn m q 
Air 2 h 3 



d 3 K ,1 



(2tt) 3 



{U 2 J q n }-Uf™U) S( 



h 2 n 2 

2m* 



2„2 



K 2 q 
2m* 



dn'{-{u 2 j q n }~uffu) = 



n M nm*q ,1 
irh 3 [ 2 



{U 2 J q n }-Uf q "U) 



^(^{r 2 ,/; i }-r/^r) 

n Mn \U\ 2 m*q 



(A4) 



irh 3 



and U =\U\T 



We have used the notation f q n = f m (q,Q'). \U\ 2 is a scalar and T is a dimensionless matrix, which is written as 
r = 1 + r s , with r s = 1/2 r s • er. Notice that T s has angular dependence because we are in the basis of eigenstates 
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of the Luttinger Hamiltonian. Thus these two contributions to the scattering term can be rewritten as 



1 rl 



n. 



Jb(/) + J m (J) = - [^{(t + r s ) 2 , /} - (i + r fl )/(i + r s )j = - [/ + {(r s + -f),f} - (/ + {r s , /} + r,/r,)] 



f - f i r 2 i i 

+ - {r s + -f, /} - - {r s , /} - - r s / r s 

T T Z T T 



r 



(A5) 



We think of r as a characteristic scattering time. The 
explicit form of the potential, determined previously, is 



U 2 = N t U 2 [l + 2a (s • S) + a 2 (s ■ S) 2 ] , 



(A6) 



where a = V/U. Everything must be averaged over the 
impurity configuration as well as directions in momen- 
tum space, and then it needs to be transformed into the 
eigenstate basis and projected onto LH and HH. When 
we do that, the term linear in a above contains only 
S z , which, when projected onto LH and HH gives some- 
thing that averages to zero over angles. Moreover, the 
configuration average of T 2 gives something which, when 
restricted to the HH and LH subspaces, is proportional 
to the identity matrix, so contributes only the scalar part 
of the scattering term. 



Mf) + Jm(f) = 



f-f , r 2 



i 



i 



— /--{r s ,/}--r s /r s . 

I T T T 

(A7) 

We separate the action of J m on the scalar and spin- 
dependent parts n, and S, of the Wigner distribution 
/ = ral + S. First on n, which is written as n = n + v , 
where v is the anisotropic part 



Jo (n) + J m (n) = (1 + T 2 ) - - - T s v - - r> 

T T T 



Averaged over impurities T 2 gives 



(A8) 



T 2 = AT, \V\ 2 [s 2 ± (S 2 + S 2 ) + 4S 2 } = , 1, (A9) 

the latter identity being valid because the matrix ele- 
ments of the S 2 restricted to the HH and LH subspaces 
are proportional to the identity matrix. We also need to 
average T s Oi T s , for which we note that Ui Uj oi = — <7j 



for i j. Averaged over impurities 
1 



r s <j x r s 



(l 2 -r 2 -r 2 )a, 



l 



F s a y F s = - (r K - r x - r 2 ) a y 



(A10) 



r s <j z v s 



- (T 2 - T 2 



T 2 y )a z 



This tells us that in the term T s S T s only the average 
of S, which we shall call S, survives. Then, writing S = 
S + E 



MS) + Jm(S) 



s r 2 - l l — = — 

- + s- (S + E) - - {F s , S} - - r s S T s 

T T T T 

(All) 

and (|A11[) we see that if we ignore 



Looking at Eq 
the term linear in r s in each of them then they do not 
mix the scalar and spin distributions. We shall work for 
now in this approximation, which is justified because the 
terms omitted are higher order in the disorder potential. 
Then we can write 



j (n) + J m (n) 



v 



1 



r 2 



T 7 T 

J_ = 1 + F 2 

T 7 



MS) + Jm(S) 



T 

e r? 



(A12) 



(S + E)--T S ST S 

T T 



The contribution to the scattering term due to the ex- 
change splitting of the bands is 



Mf) 



rq 2 h 2 



(/ - /) H pd 



m* df 
rqh 2 dq 



H 



pd 



rq 2 h 2 



ill 



(/' 



df 



(A13) 



